Monte Carlo Simulations of Conformal Theory Predictions for 
the 3-state Potts and Ising Models 



G.T. Barkema 
Institute for Advanced Study, Princeton, NJ 08540 

John McCabe* 
Laboratoire de Physique Thiorique ENSLAPP t 
ENS Lyon, 46 alii d'ltalie, 69364 Lyon Cedex 07, France 

Abstract 

The critical properties of the 2D Ising and 3-state Potts models are investi- 
gated using Monte Carlo simulations. Special interest is given to measurement 
of 3-point correlation functions and associated universal objects, i.e. struc- 
ture constants. The results agree well with predictions coming from conformal 
field theory confirming, for these examples, the correctness of the Coulomb 
gas formalism and the bootstrap method. 

I. INTRODUCTION 

Conformal field theory has produced many precise predictions for 2-dimensional, 
equilibrium, critical systems. They fall into two large groups: critical exponents |2| and 
operator-product structure constants |§. Theoretical calculations of these quantities are 
based on very special properties of the representations of the conformal group that are be- 
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lieved to be relevant to 2D critical statistical systems || - degenerate Verma modules 
and modular invariance ||. The calculation of structure constants exploits indirect tech- 
niques of constructing correlation functions - the bootstrap equation |J and the Coulomb 
gas formalism ||. For many important systems, like the Potts model, the only method avail- 
able to calculate these universal constants relies on conformal field theory. Though there is 
widespread sentiment that these techniques are correct, the abstractness of their analyses 
and the absence of independent theoretical confirmation of their predictions justify an effort 
to obtain experimental confirmation. 

Experimental tests of the predictions for the structure constants are very difficult and do 
not presently exist. They would necessitate the measurement of both 2-point and 3-point 
correlation functions at the critical point. Such measurements are much simpler in Monte 
Carlo simulations. The purpose of this article is to report on Monte Carlo experiments 
for two well-known models - the Ising and 3-state Potts models. The results will give 
"measured" structure constants that will be compared to predictions of conformal field 
theory |7|||. They will provide us with both an "experimental" test of conformal field theory 
and an insight into the methods necessary to measure these new "universal" quantities that 
it predicts in 2D critical systems. 

We will briefly summarize some facts about the Potts and Ising models that will be 
important to our analysis. A review of the statistical properties of these models can be 
found in Refs. PJlQl. Their identifications with conformal field theories are discussed in 
Refs. 0,|5],|7]|| . Since the theoretical tools necessary to perform the simulations are minimal, 
we refer to Refs. |5|,[8| for explanations of the theoretical calculation of structure constants. 

Our simulations will be for models on square lattices with periodic boundary conditions 
generated by two primitive vectors, n. The Hamiltonian of these models has the following 
form: 

H = E{x) = -J/2 [S(x + n)+ S{x - ft)} S*(x) (1) 

x x,n 

where x is a vector on the lattice. The spin density, S(x), and the energy density, E(x), are 

2 



operators that describe the coupling of the physical system to magnetic and temperature 
perturbations respectively. They are fundamental conformal operators (primary fields) hav- 
ing simple local definitions on the lattice. The field S(x) takes the discrete values ±1 for 
the Ising model and 1, exp(27ri/3), exp(— 2m /3) for the 3-state Potts model. The energy 
density operator is defined locally by the value of S(x) on five neighboring lattice sites, i.e. 

E(x) = -J^2-[S(x + n)+S(x-n)]S*(x). (2) 

n 

Both of these operators exhibit scaling behavior at the critical point. 

Traditionally, the Hamiltonian for the Potts model is written in a slightly different form, 

i.e. 

H = -JJ2$(S(x),S(x + n)). (3) 

x.it 

The Hamiltonians (0) and (|3|) are equivalent, except for an overall additive constant of — J/2 
per bond, and a scaling factor of 2/3: a pair of aligned neighboring spins contributes 1 to 
the summation in both cases, but a pair of non-aligned neighboring spins contributes -1/2 
in flip and in (g). 

The scaling behavior of 2-point correlations at a critical point defines the conformal 
dimension, A,, of a scaling field 4>i(x). For spinless fields like S(x) and E(x), it is given by: 

= T0- (4) 

True scaling fields (ex. conformal primary fields) have vanishing statistical averages at a 
critical point, || e.g. {4>i(x)) = 0. To obtain such fields, one must subtract the thermal 
averages from lattice fields with non-zero averages, like E(x). Only the subtracted operators 
obey the scaling law of (|J). The subtraction constants are not universal and are not described 
by conformal theory. This subtraction procedure must be explicitly done in any simulation. 

Finally, we mention that (|j) and all other equations for critical correlation functions 
manifestly respect the discrete symmetries of (1). The spin field, S(x), transforms under 
the discrete symmetry Z2 for the Ising and under Z3 for the 3-state Potts model. Its 
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correlation functions will obey superselection rules, at the critical point, associated with 
these symmetries. 

The conformal dimensions of S(x) and E(x) have been known for the Ising and Potts 
models for sometime [|],[5],[8],|9]] . Their explicit values are given in Table I. 

Our main interest concerns the predictions from conformal field theory for the 3-point 
correlation functions. It is well-known that the 3-point correlations of conformally invariant 
theories have the following special form Q: 

The quantities Cij k are the structure constants. Much of the revival of interest in conformal 
theories during the 1980's was associated with the realization that, in 2D, the Cafe's were new 
universal quantities different from critical exponents. More importantly, they were shown 
to be calculable from symmetry considerations alone ||. The iVj's define the normalizations 
of the 2-point functions. They are not universal and must be measured in our simulations 
of 2-point correlations before extracting the universal constants, Cij k , from (|5[). 

The calculation of the structure constants has been achieved for a large variety of minimal 
conformal models by using their special mathematical properties - the existence of null 
vectors fHH. These models are believed to describe the critical behavior of many of the 
important statistical systems. The critical point of the Ising model has been identified 
with the A 3 conformal minimal model ||. The critical point of the 3-state Potts model 
has been identified with a Z 3 symmetric version of the D 5 conformal minimal model |H§. 
The values of the structure constants resulting from theoretical calculations based on these 
identifications are summarized in table 0; the detailed calculations are found in Refs. J?H§. 

The structure constants not shown in table [Tl] vanish due to the discrete symmetries of 
the two models. The vanishing of Ceee is, however, less trivial. It results from a well-known 
fusion rule ($21 x $21 = 1) whose discovery stimulated the revival of interest in 2D conformal 
symmetry || in 1984. We will compare this prediction of conformal field theory directly to 
simulations. The fusion rules only state when structure constants are non-zero. The actual 



values for non-zero C^'s of Table |I| have been calculated by using two other tools of 2D 
conformal theory - the screened Coulomb gas formalism || and the so called "bootstrap" 
equations ||. Thus, the values of the non-zero Cijk s are a second fundamental prediction 
of conformal field theory. Experimental confirmation of their explicit values supports the 
validity of these latter two tools of 2D conformal theory. 

While the constants of the Ising model were known before the arrival of conformal field 
theories fil l, those for the 3-state Potts model have not been found by other methods |§. 



Thus, the later model allows a real test of the methods of conformal field theory. 

II. ANALYSIS 

Our simulations utilize the following procedure. First, the infinite lattice critical tem- 
perature is found from duality considerations T c = Tp. Next, the exact value of T c , for our 
finite lattice, is determined by calculating the 2-point correlations of S(x) near T^: below T c , 
the correlations approach a constant value at large distances, above T c , the correlations fall 
off exponentially, and exactly at T c the correlations show power-law scaling behavior. After 
constructing the 2-point correlation of S(x) at T c , we can measure the scaling dimension of 
S(x) and the normalization constant N$ in (§]). Next, the thermal average of E(x), (E), 
is measured at the critical temperature. Then, the critical 2-point correlation of E(x) is 
simulated. The nonscaling contribution (E) 2 is subtracted, and we then determine the scal- 
ing dimension of E(x) and Ne- Finally, the 3-point correlations are simulated. The scaling 
exponents can be extracted and the structure constants are found with the help of (|5D and 
the values of N$ and Ne- The last step is to compare our "simulated" scaling dimensions 
and structure constants with the conformal field theory predictions of Tables | and |I[ 

The Monte Carlo simulations of the Ising model were carried out on a 512 x 512 square 
lattice with periodic boundary conditions. The algorithm used to generate sample configu- 



rations is a cluster algorithm, as outlined by Wolff [12]. In one such cluster move, the time 



scale was incremented by the fraction of spins included in the cluster. In the work presented 
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here, this time scale is only relevant for defining the thermalization time and time between 
samples. In our Ising simulations, we thermalized over 300 such time units and took 18000 
samples separated by 20 time units. To obtain statistical error bars, these samples were 
blocked in groups of 150 samples, and their standard deviation was obtained. 

By determining the temperature at which the 2-spin correlation (S(x)S(O)) shows power- 
law scaling, we found that for our system size T c /Td = 1.0015 ± 0.0002. All data reported 
in this work for the Ising model is obtained with this lattice size and at the above value 
for T c . There are two independent theoretical calculations of Cses for the Ising model, and 
one does not rely on conformal theory [|TTJ. Here, we are sure of the value of Cses and the 



Ising model is a test of our procedure. It gives us more assurance when simulating the Potts 
model where no independent checks exist. The simulations of the Potts model are the real 
experimental tests of conformal field theory. 

Fig. HI shows the 2-point correlations for the Ising model. It is of particular importance 
to mention that the vacuum expectation value, (E) 2 , has already been subtracted in graphs 
of correlation functions of E(x). The measured exponents rjss = 0.26 ± 0.01 and t)ee — 
2.01±0.05 are in agreement with the values presented in Table |, i.e. four times the conformal 
dimensions for each field. More importantly, Fig. [I] tells us that iVj = 0.704 ± 0.004 and 
N 2 E = 0.42 ±0.02. 

The 3-point correlations were measured by placing one operator at = (0, 0) (the center 
of the lattice), and the two others at f\ = (r, 0) and r 2 = (0,r), i.e. two vectors along the 
two perpendicular lattice directions at a distance r. Fig. [| shows this correlation for the 
Ising model as a function of r. Conformal field theory, eq. ([5]) and Table |, predicts that for 
the Ising model: 

(S n E Q S n ) = 2^N 2 s N E C S Esr- b/i . (6) 

The observed power law behavior t]ses — 1-30 ± 0.05 is in agreement with the predicted 
exponent. The prefactor ksES = 0.33 ±0.02 and the measured values of Ns and Ne allow us 
to calculate the structure constant. We find that Cses = 0.54 ± 0.05. This measured value 



compares well with the theoretical value of Table |T[ Having shown that our simulations work 
for the Ising model, we can proceed to simulating the Potts model with some confidence. 

For the 3-state Potts model, we repeat the analysis. The Monte Carlo simulations were 
performed on a 500 x 500 lattice with the Wolff algorithm. Again, we thermalized over 300 
time units, and took samples separated by 20 time units. For the Potts model, the total 
number of samples was 21000. As for the Ising model, these samples were blocked in groups 
of 150 to obtain statistical errors. 

— # 

By determining the temperature at which the 2-spin correlation (S(x)S*{0)) shows 
power-law scaling, we found that for our system size T c /To = 1.0005 ± 0.0003. The sim- 
ulations of the 2-point functions give JVj = 0.54 ± 0.03 and N% = 0.125 ± 0.005 (Fig. g). 
The power-law dependence of these correlations were measured as r]ss* = 0.26 ± 0.02, in 
agreement with Table |, and (rjss* = 1/4), and t\ee = 1-66 ± 0.04, in slight disagreement 
with Table | (tjee = 8/5). If we again place one operator at = (0,0) (the center of the 
lattice) , and the two others at r\ = (r, 0) and r 2 = (0, r) , the values of Table | in combination 
with (|5|) predict that for the Potts model: 

(S ri S S r2 ) = 2- 1 ^NlC SSS r- 2/5 . (7) 

and 

(S ri E S; 2 ) = 2^NlN E C S ES*r- 1&/l \ (8) 

The measurements shown in Fig. f| fit this power-law behavior well: (S ri SoS r2 ) = 
k S ssr~ vsss , where k SS s = 0.44±0.04 andr/ 55S = 0.39±0.02, and (S ri E SZ 2 ) = k SE s*r- VsES * , 
where Icses* = 0.14 ± 0.01 and t/ses* = 1 11 ± 0.04, respectively. The exponents agree 
with the theoretical values 2/5 and 16/15. Combining these results with fl5|), we ob- 
tain measured values for the two non-zero structure constants of the 3-state Potts model, 
Cses* — 0.61 ± 0.06 and Csss — 1.16 ± 0.14. Again, the agreement with the predictions of 
Table [TJ is quite good. 

Finally, we show in Figs. [|, |5| the simulations of the 3-point correlation (E ri EoE r2 ) . For 
both the Ising and the Potts models, this correlation should vanish due to the conformal 
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fusion rule $21 x $21 = 1 @H- As these graphs show, the 3-point correlation functions 
are consistent with zero when the distances between the operators are as small as 5 lattice 
spacings (Note that the scale of the graphs is greatly magnified.). At very short distances 
the correlations do not vanish, because the system sees the discreteness of the lattice and 
is thus not described by conformal theory. The vanishing of these correlations at distances 
greater than a few lattice spacings is strong evidence for the correctness of the conformal 
fusion rule. 

In summary, we have measured 2- and 3-point correlations for the two-dimensional Ising 
and 3-state Potts models and have compared both their exponents and prefactors, e.g. 
structure constants, with predictions from conformal field theory. All measured exponents, 
except one, are within one standard deviation of theoretical predictions, and the remaining 
one is within two standard deviations. Far more interesting is the fact that the structure 
constants are also in good agreement with conformal theory predictions. For Ceee, this 
gives us a verification of one fusion rule of conformal models. For Cses* an d Csss, this gives 
a direct test of the screened Coulomb gas and bootstrap equation formalisms which were 
necessary to obtain their theoretical values. Our Monte Carlo simulations strongly support 
the validity of the detailed conformal theory methods that have allowed the calculation of 
higher correlation functions at critical points. 
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Monell Foundation. JM wishes to thank T. Wydro for helpful discussions and to also thank 
the Laboratoire de Physique Theorique of the ENS de Lyon for its hospitality during initial 
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TABLES 
TABLE I. Scaling dimensions of 

the Spin and Energy Density Fields 



Model: 


Ising 


3-state Potts 


S(x) 


1/16 


1/15 


E(x) 


1/2 


2/5 



TABLE II. Predictions of structure constants 
from conformal theory 

Model Structure constant value 



Ising C S es 1/2 

Ceee 

3-state Potts C S es* 0.546 

Ceee 



Csss (= Cs*s*s*) 1-092 
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FIGURES 

FIG. 1. 2-spin correlations (top) and energy density 2-point correlations (bottom) in the Ising 
model on a 512 x 512 lattice at T = T c . The dashed lines are our fit to the power-law behavior, given 
by: (S S r ) = k S sr~ vss , where k SS = 0.704±0.004 and r] SS = 0.26±0.01, and (E E r ) = kEE^~^ EE i 
where k E E = 0.42 ± 0.02 and r] EE = 2.01 ± 0.05. 

FIG. 2. Mixed 3-point correlations of the Ising model on a 512 x 512 lattice at T = T c . The 
dashed line is our fit to the power-law behavior, given by: (S r2 E S r2 ) = ksES r ~ VSES > where 
k S ES = 0-33 ± 0.02 and r] SES = 1.30 ± 0.05. 

FIG. 3. 2-spin correlations (top) and energy density 2-point correlations (bottom) in the 
3-states Potts model on a 500 x 500 lattice at T = T c . The dashed lines are our fit to the power-law 
behavior, given by: (S S*) = k S s* r'Vss* ) where k S s* = 0.54 ± 0.03 and rjss* = 0.26 ± 0.02, and 
(E E r ) = k E Er~ VEE , where k E E = 0.125 ± 0.005 and t] E e = 1-66 ± 0.04. 

FIG. 4. 3-spin correlations (top) and mixed 3-point correlations (bottom) in the Potts model 
on a 500 x 500 lattice at T = T c . The dashed lines are our fit to the power-law behavior, 
given by: (S ri S S r2 ) = k S ssr~ Vsss , where k SS s = 0.44 ± 0.04 and r] SS s = 0.39 ± 0.02, and 
(S ri E S* 2 ) = k S ES*r-vsEs* j w here k SE s* = 0.14 ± 0.01 and 7] S es* = 1-H ± 0.04. 

FIG. 5. Energy-density 3-point correlations of the Ising model on a 512 x 512 lattice at T = T c . 
For r > 10, (E ri E E r2 ) < 2 • 10~ 4 . 

FIG. 6. Energy-density 3-point correlations of the Potts model on a 500 x 500 lattice at T = T c . 
For r > 10, (E ri E E r2 ) < 10~ 5 . 
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